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Abstract 

We use the shell model Monte Carlo approach to study thermal properties 
and pair correlations in 54,56,58pg ^.^^ ^^Cr. The calculations are performed 
with the modified Kuo-Brown interaction in the complete IpO/ model space. 
We find generally that the proton-proton and neutron-neutron J = pairing 
correlations, which dominate the ground state properties of even-even nuclei, 
vanish at temperatures around 1 MeV. This pairing phase transition is accom- 
panied by a rapid increase in the moment of inertia and a partial unquench- 
ing of the Ml strength. We find that the Ml strength totally unquenches 
at higher temperatures, related to the vanishing of isoscalar proton-neutron 
correlations, which persist to higher temperatures than the pairing between 
like nucleons. The Gamow- Teller strength is also correlated to the isoscalar 
proton-neutron pairing and hence also unquenches at a temperature larger 
than that of the pairing phase transition. 

PACS numbers: 21.60.Cs, 21.60.Ka, 27.40.-hz, 21.10.Ma 



1 



Typeset using REVTgX 



I. MOTIVATION 



It has long being recognized that pairing correlations play an essential role in low-energy 
nuclear physics. In fact, mean-field theories IQ like Hartree-Fock+BCS (which adds pairing 
as an additional degree of freedom to the mean-field solution) or the more consistent Hartree- 
Fock-Bogoliubov theory have proven to be quite successful in the general description of 
ground states of even-even nuclei. Although in many cases pairing among like nucleons 
is sufficiently well approximated by Cooper pairs (in which a nucleon is coupled to its 
time-reversed partner) phenomenology indicates that d-wave pairs of like nucleons play an 
important role in deformed heavier nuclei, such as those of the rare-earth region. This 
observation prompted the development of the Interacting Boson Model, which is remarkably 
successful in describing low-lying nuclear spectra 0. 

The properties of nuclei at finite temperature have attracted recent experimental 
and theoretical interest. Although it is obvious that the influence of pairing decreases 
with temperature and must vanish in the high-temperature limit, the mean- field description 
of nuclei at finite temperature is inadequate due to the neglect of important quantum and 
statistical fluctuations 0. Recently developed shell model Monte Carlo (SMMC) methods ^ 
do not have this shortcoming and also allow the consideration of model spaces large enough 
to account for the relevant nucleon-nucleon correlations at low and moderate temperatures 
(T < 1.4 MeV [|). 

A first study considered the thermal properties of ^^Fe within the SMMC approach . 
An important result was that the like-pair correlations, described by the BCS-like proton 
and neutron monopole pairs, vanish in a small temperature interval around T = 1 MeV. 
This pairing phase transition is accompanied by a rather sharp rise of the moment of inertia 
and a partial unquenching of the B[M1) strength. However, we also found that the total 
Gamow- Teller strength B[GT^) remains roughly constant through the phase transition and 
does not unquench. It was conjectured [§ that the isoscalar proton-neutron correlations, 
which are known to be the main source for the quenching of the B{GT+) strength P,p!0|, 
persist to higher temperatures than the like-pair correlations. 

The purpose of this paper is to demonstrate the validity of that conjecture. To that 
end we define improved measures for the strength of the various pair correlations and study 
their temperature dependence, as well as their correlation with such nuclear observables 
as the moment of inertia, and the total B{M1) and B{GT+) strengths. Detailed SMMC 
calculations are presented for the nuclei 54,56,58pg ^^Cr. 

Our paper is organized as follows. After a brief discussion of the SMMC method in 
Section 2, we introduce our definition of the pairing correlation strengths. In Section 3 we 
present the results of our study. We begin by discussing the temperature dependence of the 
moment of inertia, the B{M1) and B{GT+) strengths, and the monopole pairing fields for all 
four nuclei, demonstrating that the results found in Ref. |^ appear to be a general feature 
of even-even nuclei with A ^ 56. In Section 3.B we study the temperature dependence 
of the various pair correlations in the four nuclei, paying particular attention to the phase 
transition related to the breaking of monopole pairs as found in [^] and the presence of 
proton-neutron correlations at higher temperatures. Finally we relate in Section 3.C the 
calculated temperature dependences of the moment of inertia and the B{M1) and B{GT^) 
strengths to the temperature dependences of selected pairing correlations. 
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II. THEORETICAL BACKGROUND 



We exploit recently developed Monte Carlo techniques |7|JTT| to calculate the thermal 
properties of 54,56,58pg 56 Qj, ^ complete Qhu model space with a realistic interaction. 
The methods we use describe the nucleus by a canonical ensemble at temperature T = (3~^ 
and employ a Hubbard-Stratonovich linearization of the imaginary-time many-body propa- 
gator, e~^^ , to express observables as path integrals of one-body propagators in fluctuating 
auxiliary fields 0,0. Since Monte Carlo techniques avoid an explicit enumeration of the 
many-body states, they can be used in model spaces far larger than those accessible to 
conventional methods. The Monte Carlo results are in principle exact and are in practice 
subject only to controllable sampling and discretization errors. To circumvent the notori- 
ous "sign problem" encountered in the Monte Carlo shell model calculations with realistic 
interactions |T2|, we adopted the procedure suggested in Ref. ^ which is based on an ex- 
trapolation from a family of Hamiltonians that is free of the sign problem to the physical 
Hamiltonian. 

Our calculations include the complete set of 1^3/2,1/20/7/2,5/2 states interacting through 
the realistic Kuo-Brown Hamiltonian [^] as modified in |T^. Some 5 x 10^ configurations of 
the valence neutrons and protons moving in these 20 orbitals are involved in the canonical 
ensemble. As in Ref. the results presented below have been obtained in MC shell model 
studies with a time step of A/? = 1/32 MeV~^ using about 5000 independent Monte Carlo 
samples at six values of the coupling constant g spaced between —1 and and the value 
X = 4 in the decomposition of the Hamiltonian. A linear extrapolation to the physical case 
{g = 1) has been adopted for the observables discussed below. 

The results presented in this paper correspond to various nuclear observables at selected 
temperatures: the total 5 (Ml) and Gamow- Teller strengths and the moment of inertia 
/, given by / = 3{.P)/T. The total B{M1) strength is defined as B{Ml) = {fl'^), where 
the magnetic moment /x is given by /2 = J2i I-I'N \^9d + gsS^, A* at is the nuclear magneton 
and gi,gs are the free gyromagnetic ratios for angular momentum and spin, respectively 
{gi = l,gs = 5.586 for protons and gi = 0,gs = —3.826 for neutrons). The total Gamow- 
Teller strength is given by B{GT^) = {{ar+y). We also explore the isovector monopole 
pairing, as described in Section 3. A. 

The main focus of this paper is on pairing correlations and their relation to nuclear 
observables. In our complete Ohuj /p-shell model space, a pair of protons or neutrons with 
angular momentum quantum numbers [JM) is generated by (a = vr for protons and a = u 
for neutrons) 



^JAf (Ja, Jb) = , \ [at 41 , ,,,, (1) 



1 + 



Jb 



(JM) 



where vrj (i/j) creates a proton (neutron) in an orbital with total spin j. Proton-neutron pairs 
can have either isospin T = 1 or T = 0. The respective definitions for the pair operators are 

AjAlUa, Jb) = , ^ , + -r^jA] (2) 



2^1 + 

for the isovector proton-neutron pairs and 
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AnUa, Jb) = ^ I , ^ - AA] , (3) 



3a3h 



(JM) 



for the isoscalar pair. 

A general description of a pair, A^jj^^, is then given by a superposition of the pair opera- 
tors: 

^JM= f^JM{ja,jb)AlM{ja,jb)- (4) 

ja>jb 

In mean-field theories like Hartree-Fock-|-BCS or HFB, the weights /3 can be related to the 
two-body potential matrix elements by the variational minimization of the free energy in 
the mean-field model space. However, such a relation is not to our disposal in the complete 
shell model and thus the weights are apriori undetermined. To overcome this arbitrariness 
of the definition (4) we will define the pairing strength as follows. With (1-3) we build a 
generalized pair matrix 

ML' = E(4M(ja, Jb)AjM{jc: 3d)) (5) 

M 

which corresponds to the calculation of the canonical ensemble average of two-body operators 
like 7r|7r27r37r4. The index a distinguishes the various possible {ja,jb) combinations (with 
ja > jb)- For example, the square matrix M has dimension A^j = 4 for J = 0, A^j = 7 for 
J = 1, A^j = 8 for J = 2, 3. In the second step, the matrix M"^ is diagonalized (/3 = 1, Nj), 
labelling the eigenvectors and eigenvalues in decreasing order 

M'xfs = (6) 

The presence of a pair condensate in a correlated ground state will be signaled by the largest 
eigenvalue for a given J, Xf, being much greater than any of the others. 

To discuss the temperature dependence of pairing correlations it is convenient to intro- 
duce an overall measure for the pairing strength. We define the pairing strength for pairs 
with spin J as the sum of the eigenvalues of the matrix M-^, 

^' = EA^ = E^ia- (7) 

(3 a 

We note that for a model space spanned by a single j-shcU one finds a simple sum rule for 
proton and neutron pairs 

^P^ = iAr„(Ar,-i), (8) 
J ^ 

where A^^, is the number of valence protons or neutrons in the shell. This sum rule holds 
approximately at low temperatures to the protons of the nuclei studied in this paper, as the 
valence protons then occupy mainly the jjj2 orbital. The sum rule is not valid at higher 
temperatures as the valence particles are then spread out over several orbitals. 

With our definition (7) the pairing strength is non-negative and indeed positive at the 
mean-field level. As we will focus on the pairing correlations beyond the mean field in this 
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paper, we also calculated the mean-field pairing strength, P^^ by the same procedure as 
outlined above, however, replacing the expectation values of the two-body matrix elements 
in the definition of the generalized pair matrix M"^ by 

(440304) nin2 - fe^u) (9) 

where = (a^a^.) is the occupation number of the orbital k. 

We then finally define the pairing correlation in the nucleus as the pairing strength 
beyond the mean field, 

-^corr = ~ Pmi- (lO) 



III. RESULTS 

A. Signals for the phase transition 

Ref. [H reported on SMMC calculations of the thermal properties of ^'^Fe within the 
complete p/-shell. One quantity studied was the thermal response of BCS-like monopole- 
pair correlations among protons and neutrons, defined as (a = vr, i/ for protons and neutrons, 
respectively) 

^BCS = ^m^ln- (H) 

m>0 

It was observed that these BCS-like pairs break at temperatures near 1 MeV. Of the ob- 
servables studied in Ref. |^ three exhibit a particularly interesting behavior with increasing 
temperature near the phase transition: a) the moment of inertia / rises sharply; b) the Ml 
strength B{M1) shows a sharp rise, but unquenches only partially; and c) the Gamow- Teller 
strength B{GT^) remains roughly constant and strongly quenched. 

We will now show that the breaking of BCS-like pairs at temperatures near 1 MeV and 
the related behavior of the three observables (moment of inertia and Ml and Gamow- Teller 
strengths) is not a particular result for ^^Fe, but is a more general feature of even-even nuclei 
in the A = 56 mass region. To this end, we present the results of SMMC calculations for 
^^Cr and 54,56,58pg selected temperatures. Fig. 1 shows the temperature dependence of 
/, B{M1) and B{GT^) for all four nuclei. Additionally, the expectation values for BCS- 
like pairing fields, (Ag^^gABcs)) for both protons and neutrons are plotted as functions of 
temperature. It is obvious from Fig. 1 that the like-pair correlations disappear at around 
1 MeV for all four nuclei. This phase transition is accompanied by the same behavior of 
/, B{M1) and B{GT^) as discussed in Ref. [H and summarized in points a)-c) above. We 
note that Ref. employed a different residual interaction (the Brown-Richter interaction 
||T5| ) than used here; the temperature dependence of the quantities shown in Fig. 1 are 
apparently largely independent of the interaction. 

While Ref. § presented results only up to T = 2 MeV, we have continued the calculations 
in this paper to higher temperatures. Although the quantitative results are expected to be 
affected by our finite model space at temperatures larger than about 1.4 MeV 0, the 
qualitative features of the observables are likely still meaningful. The moment of inertia 
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shows a T~^-decrease with temperature at T > 2 MeV, as is expected from the definition 
/ = 3(J^)/T; our finite model space requires (J^) to approach a constant in the high- 
temperature hmit. The B{M1) strength unquenches in two steps. Following its partial 
unquenching at the pairing phase transition, it remains roughly constant up to T = 2.6 
MeV, where it finally starts to fully unquench. Although not shown in Fig. 1, we have 
checked that in the high-temperature limit the B{M1) strengths approach the appropriate 
values. Being roughly constant at lower temperatures, B{GT+) unquenches at T > 2.6 MeV. 
We observe that this happens simultaneously with the second step of the unquenching of 
B{M1) and, as we will see below, the two phenomena have a common origin. Again we have 
checked that B[GTj^) approaches the appropriate values in the high-temperature limit. 



B. Pair correlations 

We have studied the pair correlations in the four nuclei for the various isovector and 
isoscalar pairs up to J = 4. Detailed calculations have been performed for selected temper- 



atures between T = 0.5 MeV and 8 MeV. As has been shown in Refs. [0,^, the SMMC 
studies at T = 0.5 MeV approximate the ground state properties for even-even nuclei well. 
We have checked that in the high-temperature limit, which in our calculations corresponds 
to T = 1/ (A/3) = 32 MeV, the pairing correlations approach the appropriate mean-field 
values. 

Although it is the pairing correlation P^orr (i-^-^ pairing strength relative to the mean- 
field) that matters for the behavior of the physical observables, we will at first discuss some 
features of the pair matrix M and its spectrum. In Fig. 2 we show the pair matrix eigenvalues 

for the three isovector J = and the isoscalar l"*" pairing channels as calculated for the 
iron isotopes ^^~^^Fe. We compare the SMMC results with those derived on the mean-field 
level, as discussed in section 2. Additionally, Fig. 2 shows the diagonal matrix elements of 
the pair matrix Ma^a, where we use the notation a = 1, ..,4 for J = 0^ pairs in the /7/2, 
P3/2, /5/2 and pi/2 orbitals. For the isoscalar pairs only the three largest diagonal matrix 
elements of M are shown corresponding to pairs in (/7/2)^ (« = 1), (^3/2)^ (« = 2) and 
if 7/2/5/2) (« = 3) proton-neutron configurations. As expected, the protons occupy mainly 
/7/2 orbitals in these nuclei. Correspondingly, the {A^A) expectation value, Mn, is large 
for this orbital; the other diagonal elements of M are small. For neutrons, the pair matrix 
is also largest for the /7/2 orbital. The excess neutrons in ^^.sspg Qgc^py i^j^e p^^2 orbital, 
signalled by a strong increase of the corresponding pair matrix element M22 in comparison 
to its value for ^^Fe. Upon closer inspection we find that the proton pair matrix elements 
are not constant within the isotope chain. This behavior is mainly caused by the isoscalar 
proton-neutron pairing. The dominating role is played by the isoscalar l"*" channel, which 
couples protons and neutrons in the same orbitals and in spin-orbit partners. We thus find 
for ^^'^^Fe that the proton pair matrix in the /5/2 orbital, M33, is larger than in the P3/2 
orbital, although the latter is favored in energy. For ^^Fe, this ordering is inverted, caused 
by the increasing number of neutrons in the P3/2 orbital which in turn also increase the 
proton pairing in this orbital. 

After diagonalization the proton pairing strength is essentially found in one large eigen- 
value. Furthermore we observe that this eigenvalue is significantly larger than the largest 
eigenvalue on the mean-field level supporting the existence of a proton pair condensate in 
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the ground state of these nuclei. For the neutrons the situation is somewhat different. For 
^^Fe, httle additional coherence is found beyond the mean-field value, reflecting the closed- 
subshell neutron structure of this nucleus. For the two other isotopes, the neutron pairing 
exhibits two large eigenvalues. Although the larger one exceeds the mean-field value and sig- 
nals noticeable additional coherence across the subshells, the existence of a second coherent 
eigenvalue shows the limits of the BCS-like pairing picture. 

In Fig. 3 we compare the SMMC pairing strengths (Eq. (7)) for the ^^Fe ground 
state (calculated at T = 0.5 MeV) with the mean-field values. Although the largest values 
for are found for pairs with larger J- values (e.g. J = 2 and 4 for isovector pairs), 
these values simply reflect the larger combinatorial possibilities to make these J-pairs in our 
model space, as the mean-field values are close to the SMMC values, i.e. there are no true 
correlations. Fig. 3 also exhibits the expected odd-even staggering: in the isovector pairing 
channels, pairs with even J are more likely than those with odd J, while it is vice versa for 
isoscalar pairs. It is also obvious from Fig. 3 that the most significant physical difference 
between the SMMC and mean-field pairing occurs in the isovector J = proton-proton 
and neutron-neutron channels. Here, the SMMC calculations exhibit a significant excess of 
pairing correlations over the mean field, reflecting the well known coherence in the ground 
states of even-even nuclei. As the nucleus ^^Fe is semimagic for neutrons, the excess is larger 
for protons than for neutrons. Due to the sum rule (8), which is approximately fulfilled for 
the ^^Fe ground state, the excess of J = pairs is counterbalanced by correlation deficiencies 
in the other isovector pairs. 

In the following we will discuss the physically important pairing correlation as defined 
in Eq. (10). Fig. 4 shows the temperature dependence of the pair correlations for selected 
pairs, which, as we will see in the next subsection, play an important role in understanding 
the thermal behavior of the moment of inertia and the total Ml and Gamow- Teller strengths. 

The most interesting behavior is found in the J = proton and neutron pairs. As 
mentioned above, the large excess of this pairing at low temperatures reflects the ground 
state coherence of even-even nuclei. With increasing temperature, this excess diminishes, 
vanishing at around T = 1.2 MeV. This behavior is in agreement with the pairing phase 
transition deduced in Ref. 0] and above from the temperature dependence of the BCS- 
like monopole pairs. We observe further from Fig. 4 that the temperature dependence 
of the J = proton-pair correlations is remarkably independent of the nucleus, while the 
neutron pair correlations show interesting differences. First, the correlation excess is smaller 
in the semimagic nucleus ^^Fe than in the others. Among the iron isotopes, the neutron 
J = correlations vanish at higher temperatures with increasing neutron number. More 
quantitatively, we have fitted the correlation excess to a Fermi function. 



We then find (in MeV) Tq = 0.82 ± 0.13, AT = 0.13 ± 0.07 for ^^Fe, Tq = 0.96 ± 0.03, AT = 
0.11 ± 0.02 for s^Fe, and Tq = 1.07 ± 0.04, AT = 0.18 ± 0.04 for ^^Fe, while the fit for ^^Cr 
is very similar to that for ^^Fe, Tq = 1.06 ± 0.04, AT = 0.14 ± 0.03. Similar fits to the 
J = proton pairing excess results in Tq ^ 0.9 ± 0.04, AT = 0.12 ± 0.03 for all nuclei. The 
vanishing of the J = proton and neutron pair correlations is accompanied by a significant 
increase in the correlations of the other pairs. 




(12) 
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The isovector J = proton-neutron correlations arc positive at all temperatures for all 
four nuclei, with a slight excess at low temperatures that increases by about a factor 3 after 
the J = proton and neutron pairs have vanished. The correlation peak is reached at higher 
temperatures with increasing neutron number (about 1.3 MeV for ^^Fe, 1.6 MeV in ^^Fe), 
while the peak height decreases with neutron excess. 

The isoscalar proton-neutron J = 1 pairs show an interesting temperature dependence. 
At low temperatures, when the nucleus is still dominated by the J = proton and neutron 
pairs, the isoscalar proton-neutron correlations show a noticeable excess. But more inter- 
estingly, they are roughly constant and do not directly reflect the vanishing of the J — 
proton and neutron pairs. However, at T > 1 MeV, when the proton and neutron-pairs are 
broken, the isoscalar J = 1 pair correlations significantly increase and have their maximum 
at around 2 MeV, with peak values of about twice the correlation excess in the ground 
state. A similar (but milder) increase is observed in the other isoscalar correlations. In 
contrast to the isovector J — proton-neutron pairs, the correlation peaks are reached at 
lower temperatures with increasing neutron excess. We also observe that these correlations 
die out rather slowly with temperature. 

In summary, our SMMC calculations of nuclei around A = 56 support the following 
temperature hierarchy of pairing correlations. In the ground state and at low temperatures, 
J — proton and neutron pair correlations dominate. When these vanish at temperatures 
around 1 MeV, proton-neutron correlations become more important, with the isovector 
proton-neutron correlations, in our model space, vanishing at temperatures around 2 MeV. 
At high temperatures, isoscalar proton-neutron correlations dominate. 



In this subsection we will discuss the relationship between pairing correlations and the 
temperature response of the Ml, and Gamow- Teller strengths, as well as of the moment of 
inertia. As in our definition of pairing, we will focus on the behavior of these observables 
beyond the mean-field values. In terms of the neutron and proton occupation numbers n„ 
and Up, respectively, the mean- field values for these observables are given by 



C. Pairing correlations and observables 




(13) 



B{M1U = E [\mp)\j)Wij)npii) + \mn)\j)\'nnij)nnii) 



(14) 




(15) 



where the blocking factors are 





(16) 
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and fl{p) and fl{n) are the proton and neutron magnetic moments. 

The SMMC results for all three observables are significantly smaller than the mean-field 
values. To quantify this suppression, we introduce quenching strengths 

(-^) quench = (-^)mf " {I)', (18) 

5(Ml)q„ench = 5(Ml)^f - 5(M1); (19) 
quench 

= B{GT+U-B{GT+). (20) 

We expect the quenching strength of the observables to be larger when the relevant pair- 
ing correlations are. Fig. 5 shows the temperature dependence of the quenching strengths 
for the three observables and relates it to selected pairing correlations, where the latter are 
given in arbitrary units. The following observations are evident and apply to all four nuclei: 

The quenching of the moment of inertia shows a sharp drop at around T = 1 MeV and 
then slowly diminishes; the SMMC and mean-field values are identical for T > 2 MeV. The 
drop at T ^ 1 MeV clearly signals a relation to the pairing phase transition observed in 
the J = proton and neutron correlations (see Fig. 4). From the definition / ~ (J^) = 
((J„, + JpY), where Jp, J„ are the proton and neutron angular momenta, respectively, we 
expect that the moment of inertia is sensitive to the proton, neutron and proton-neutron 
correlations. This expectation is confirmed in Fig. 5, which shows that the temperature 
dependence of the quenching is nearly identical to the sum of the isovector J = pair 
correlations. The steep drop at T ^ 1 MeV is thus due to the J = proton and neutron pair 
correlations. We note that the moment of inertia unquenches more slowly with temperature 
than the proton and neutron J = pairs. This difference is due to the isovector J = 
proton-neutron correlations. 

The quenching of the Ml strength decreases significantly near T = 1 MeV. Following 
a minimum at T 1.3 MeV, the quenching increases again and then, after a maximum 
near 2 MeV, decreases slowly. While the drop at T = 1 MeV again signals an association 
with the J = proton and neutron pairing phase transition, the maximum at 2 MeV as 
well as the slow dying out clearly resembles the temperature dependence of the isoscalar 
J = 1 proton-neutron correlations. These two correlations reflect the two components in 
the Ml strength. The orbital part is sensitive to the J = proton pairing correlations (the 
gyromagnetic moment gi is zero for neutrons) , while the quenching of the spin component is 
dominated by isoscalar proton- neutron correlations PJTO[]. This conjecture is in fact nicely 
confirmed, as for all four nuclei the temperature dependence of the quenching of the Ml 
strength is well described by that of the sum of the J = proton and isoscalar J = 1 
proton-neutron correlations. 

Our calculation reveals a close (linear) relation between the orbital part of the Ml 
strength and the moment of inertia at low temperatures, where both quantities are domi- 
nated by isovector J = pairing correlations. Such a linear relation is experimentally known 



for heavier nuclei and is supported by QRPA calculations |1T7HT9|]. 

The quenching of the B{GT^) strength is roughly constant at low temperatures (T < 1 
MeV). It then increases slightly, before it slowly dies out following a maximum at T ^ 2.5 
MeV. Fig. 5 shows that the temperature dependence of the Gamow- Teller strength is driven 
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by the isoscalar proton-neutron correlations, as the rough constancy at low temperatures, the 
maximum at around T = 2.5 MeV and the slow decrease of the correlations towards higher 
temperatures of the B{GT^) quenching qualitatively mimics the temperature dependence 
of the isoscalar J = 1 proton-neutron correlations. With increasing neutron excess the 
maximum in the B{GT^) quenching at T ^ 2.5 MeV is weakened, in qualitative agreement 
with the J = 1 proton-neutron correlations. 

IV. SUMMARY 

In a previous SMMC study of the nucleus ^^Fe, a sharp rise in the moment of inertia and 
the Ml strength at a temperature T ^ 1 MeV were related to the vanishing of the BCS- 
like proton and neutron Cooper pairs. In this paper we have investigated in more details 
this phase transition and its connection to the temperature dependence of the moment 
of inertia and of the Ml and Gamow- Teller strength. To do so we have performed SMMC 
calculations for the iron isotopes 54,56,58pg ^^^^ ^^Cr within in the complete p/-shell using the 
KB3 interaction. To study the pairing correlations within the shell model, we have defined 
a generalized pair matrix for isovector and isoscalar pairs for a given angular momentum 
and have proposed the sum of the eigenvalues of this matrix (trace) as a measure of the 
pairing strength. Using the same formalism, but replacing the two-body pair operators in 
the generalized pair matrix by occupation numbers we have derived the mean-field pairing 
strength. The physically relevant pair correlations are then defined as the difference of the 
SMMC and mean-field pairing strengths. 

A detailed investigation of the correlations in the various pairing channels gives the 
following qualitative picture. 

At low temperatures the SMMC calculations show a strong excess of J = proton and 
neutron pairs over the mean-field values. This excess drops sharply at around T = 1 MeV 
for all four nuclei studied here, confirming the pairing phase transition conjectured in Ref. 
1^. After the vanishing of the dominating J = pairs, proton- neutron correlations show 
an excess over the mean-field values. We find that isoscalar proton-neutron correlations 
persist to higher temperatures than the isovector correlations and, within our model space, 
vanish slowly at temperatures larger than 3 MeV. We observe some evidence that the pairing 
correlations depend on the neutron excess of the nucleus. This nuclear isospin-dependence 
warrants more detailed explorations. 

The temperature hierarchy of the pairing correlations and the J = proton and neu- 
tron pairing phase transition manifests itself in the temperature dependence of the three 
observables we studied in detail. The moment of inertia of the nucleus correlates with the 
sum of the isovector J = pair correlations. In particular, we find that the pairing phase 
transition is connected with a sharp rise of the moment of inertia. Due to its orbital and 
spin components, the Ml strength unquenches generally in two steps: the orbital part un- 
quenches at the J = proton and neutron pairing phase transition, while the spin part 
slowly unquenches at higher temperatures, in concert with the decrease of the isoscalar 
proton-neutron correlations. The temperature dependence of the Gamow- Teller quenching 
B[GT+) is well approximated by that of the isoscalar J = 1 proton- neutron pairs. 

Although the present calculations have been performed in model spaces larger than any 
imagined tractable until only recently, the quantitative results suffer from finite model space 
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limitations at temperatures greater about 1.5 MeV. Although we believe that the qualitative 
results presented in this paper will persist in calculations performed in even larger model 
spaces, such work is certainly warranted for a quantitatively reliable description of the 
unquenching of the isoscalar correlations at high temperatures. Such calculations are feasible 
within the SMMC approach and, as a next step, we are planning studies in the complete 
(pf+sdg) shells. 
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FIGURES 



FIG. 1. Temperature dependence of the moment of inertia (left top), the B{M1) strength 
(left bottom), the B{GT^) strength (right top), and the expectation values of the BCS proton and 
neutron pairing strength (right bottom). SMMC results are shown for ^^Fe (part a), ^^Fe (part b), 
^^Fe (part c), and ^^Cr (part d). 

FIG. 2. Diagonal elements (upper) and eigenvalues (middle) of the pair matrix M in the three 
isovector J = 0+ and the isoscalar J = 1+ channels, as calculated for 54,56,58pg at T = 0.5 MeV. 
Eigenvalues of the pair matrix on the mean-field level (see text) are shown in the lowest panels. 

FIG. 3. Comparison of the SMMC (P-^) and mean-field (P^f) pairing strengths in the ^^Fe 
ground state for selected pairing channels. 

FIG. 4. Temperature dependence of the pair correlations P^^r ^ defined in Eq. (10) for 
selected pairing channels. 

FIG. 5. Comparison of the temperature dependence of the quenching of the physical observ- 
ables (see Eqs. (18-20)) with related pairing correlations for ^^'^^Fe (part a) and ^^Fe, ^^Cr (part 
b). Upper: moment of inertia and the sum of the isovector J = pairing correlations. Middle: the 
Ml strength and the sum (solid line) of J = proton (dashed) and isoscalar J = 1 proton-neutron 
(dotted) correlations; Lower: the B{GTj^) strength and the isoscalar J = 1 proton- neutron corre- 
lations. The pairing correlations have been scaled arbitrarily. 
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